Exploration of semen quality analyzed by casa-mot systems of brahman bulls infected with BLV and BHV-1

Enzootic bovine leukosis virus (BLV) and bovine herpesvirus 1 (BHV-1) are very important infectious agents for the livestock industry worldwide. The present study aimed to explore the association between natural exposure to BLV and BHV-1 with sperm quality analyzed by Computer-Assisted Semen Analysis (CASA) systems. Ten sexually mature Brahman bulls, with sanitary status BLV+/BHV-1+ (n = 2), BLV−/BHV-1+ (n = 6) and BLV-/BHV-1- (n = 2) were evaluated twice, 30 days apart. Results showed that sanitary status of each bull was not associated with semen quality. It was found that the quality of the semen from the second collection was better due to the interruption of sexual rest. The evidence thus revealed that a bull infected with BLV generated good-quality contaminated semen and, therefore, that it is essential to detect contaminated seminal samples to prevent the spread of BLV. A multivariate analysis showed the presence of four sperm subpopulations in Brahman bulls that differ significantly in their kinematic patterns and with respect to sanitary status (P < 0.05), indicating that infection-free and seronegative bulls present the best kinematic parameters, which improved discrimination of sperm quality according to sanitary status. Overall, the analyses indicate that the seropositive-infected bulls with BLV and BHV-1 should be excluded from beef cattle farms.


Results
In this study, blood samples showed that two bulls were seropositive and infected with BLV and BHV-1 (BLV + / BHV-1 + ), six bulls were seropositive and latently infected with BHV-1 (BLV − /BHV-1 + ), and two bulls were negative for the presence of antibodies and antigen associated with the BLV and BHV-1 (BLV − /BHV-1 − ).All bulls maintained the sanitary status throughout the investigation.Furthermore, only one seminal sample from bull 498 was positive to BLV at 0 day but it was negative at 30 days (second seminal sample).The BLV detected in blood and seminal samples of bull were characterized as genotype 5 of BLV (strain 370322-498BlSp-CR, GenBank accession number: OP718193 (https:// www.ncbi.nlm.nih.gov/ nucco re/ OP718 193; accessed on 17 July 2023).No blood and semen samples were positive for the presence of BHV-1.
The volume per ejaculate ranged from 2 to 12 mL, whereas the mean concentration ranged from 388 to 1560 M/mL and the total motility ranged from 37.7 to 71.0%.Progressive motility (PM) and kinematics variables did not show a change in trend between evaluations (0 and 30 days) or between the sanitary status of the bulls that would allow differentiation or increase the discriminant potential on the seminal quality initially characterized based on the percentage of normal morphology and viability (Tables 1 and 2, Fig. 1).
The assessment period seems to affect seminal quality and the percentage of bulls with doubtful seminal quality.During the initial evaluation, 70% (7/10) of the bulls showed at least one of the seminal parameters below the expected limit of detection for a sample with satisfactory seminal quality.At 30 days, only 40% (4/10) of the bulls showed doubtful seminal quality and the proportion of bulls with satisfactory seminal quality increased (Tables 1 and 2, Fig. 1).
Viability and percentage of abnormal sperm were the values with the best discriminant potential to determine the quality of fresh semen samples.The descriptive analysis of frequency of different types of sperm abnormalities at 30 days showed that all bulls had less than 30% abnormalities in the ejaculate.Four types of abnormalities had the highest prevalence: 100% (10/10) of bulls had sperm with proximal droplets, which ranged from 0.47 to 9.34%; 90% (9/10) of bulls had sperm with loose heads, which ranged from 0.47 to 18.93%; 90% (9/10) had Table 1.Sanitary status, means of sperm parameters and semen quality of Brahman bulls (n = 10) at the beginning of the study.PM progressive motility, TM total motility, VCL curvilinear velocity, VSL straightline velocity, VAP average path velocity, ALH amplitude of lateral head displacement, LIN linearity of forward progression, STR straightness, WOB wobble, BCF beat-cross frequency, S satisfactory, D doubtful, U unsatisfactory. 1 Seropositive-infected bull with BLV and BHV-1. 2 Seropositive-latently infected bull with BHV-1. 3 Negative bull to BLV and BHV-1.sperm with tightly coiled tail, which ranged from 0.41 to 7.14% and, finally, 80% (8/10) had distal droplets, which ranged from 0.41 to 1.90%.There were also two rare anomalies in bulls, namely double heads and bent midpiece (Table 3, Fig. 2).
The PFA rendered two principal components (PCs) for sperm kinematic variables, explaining 84.51% of the total variance (Table 4).The two PCs are similar in total variance explained.The PC1, named "linearity", was positively correlated to progressivity parameters (LIN, STR and WOB).PC2, was named "velocity", and was positively correlated to velocity parameters (VCL and VAP) and sperm head oscillation (ALH).The straight-line velocity (VSL) was included in the two PCs (Table 4).
Most of all kinematic parameters showed significant differences between different sanitary status.This means that the motility was more linear in a sanitary status BLV − /BHV-1 − than that of BLV + /BHV-1 + or BLV − /BHV-1 + bulls (Table 5).When comparing the sanitary status, VSL, VAP, LIN, STR, ALH and BCF were higher in BLV − /BHV-1 − than in BLV + /BHV-1 + and BLV − /BHV-1 + bulls.In contrast, WOB showed no differences among bulls with different sanitary status.When comparing VSL, VAP, LIN and STR, they were higher in BLV − /BHV-1 + that in BLV + /BHV-1 + bulls.Following these results, it seems that the sanitary status of a bull can affect the velocity and linearity of sperm kinematic parameters (Table 5).
Four well-defined sperm kinematic subpopulations (SP1, SP2, SP3, SP4) were identified; they showed clear differences in most kinematic parameters (Table 6).SP1 had the lowest values in all kinematic parameters, being considered a "slow and nonlinear" subpopulation.SP2 had sperm with high values for VAP and VCL, but did not have high values for VSL, LIN and STR, and was regarded as "fast, but with less straight and linear" sperm than SP4.SP3 had high values for LIN, STR and WOB, but low values in VCL, VSL, VAP, being defined as "linear and slow" subpopulation.SP4 included spermatozoa with the highest linear trajectories (LIN, STR) and high speed (VSL, VAP) and BCF, and included fast straight, and linear sperm with a high tail beat (Table 6).
The percentage of spermatozoa of each sanitary status in reference to the total number of spermatozoa varied between the subpopulations.The sanitary status BLV − /BHV-1 + (SS2) had the highest percentage of spermatozoa (54.3%), followed by BLV − /BHV-1 + (SS1) with 23.1%, and finally by BLV − /BHV-1 − (SS3) with 22.7% of total spermatozoa analyzed.There was a higher percentage of SP4 in SS1 y SS2 (9.6% and 18.8%, respectively), but in SS3 was SP3 (7.6%) (Table 7).This means that bulls that were in SS1 y SS2 presented spermatozoa with high linear trajectories and speed.There were differences between sperm SPs within each SS (P < 0.05) with VCL and VAP exhibiting the highest differences between the sub-populations (Table 8).There were also differences between SPs for VSL and STR within SS2 and SS3 (P < 0.05), but not in SS1 (P > 0.05).For WOB there were differences between SPs only in SS2 and SS3; whereas for ALH there were differences between SPs only in SS1 and SS3.The SP distribution was not the same for each sanitary status (Table 8).When SP4 was compared between sanitary status, kinematics differences (P < 0.05) were found, for the variables VAP, LIN, STR and WOB.When SP1 was compared between sanitary status, kinematics differences (P < 0.05) were found only in WOB.The differences when comparing SP2 between the sanitary status presented a behavior similar to that shown in SP4, while the behavior of the differences in SP3 was similar to SP1.Following these results, SS1 and SS2 affect the kinematic variables of bull ejaculates.

Discussion
The results of this study showed that the percentage of normal sperm morphology and sperm viability were the variables with the greatest discriminating potential for samples with satisfactory, doubtful, and unsatisfactory seminal quality.These results agree with those reported earlier [19][20][21] but contrast with other studies which found that the percentage of abnormal spermatozoa was not a good predictor of fertility level 22 .
Our results indicate that bulls presented a spermiogram, which contained < 30% abnormal sperm, as a threshold value described in other studies [23][24][25] .However, there were three abnormal sperm morphology defects (loose head, tightly coiled tail and proximal droplet) that overcame 6% of individual thresholds for each abnormality 26 .The concept of sperm with compensable or uncompensable traits allows to differentiate thresholds as the case may be 27 .For example, compensable traits such as loose heads preclude affected sperm from fertilizing the ovum and the threshold in this case can be 30% 26 .For uncompensable traits the authors suggested that threshold of such abnormalities should not surpass 20% 26,28 .At 30 days, all the bulls showed a percentage of abnormal spermatozoa lower than 30%, a value associated with satisfactory quality for the seminal sample 29 .The prevalence of morphological abnormalities provides very valuable information to assess reproductive health associated with spermatogenesis 26 .When considering the prevalence of each sperm abnormality in each bull, the presence of proximal droplets was the most prevalent sperm abnormality, followed by loose heads, tightly coiled tails and distal droplets, which are within the 3 major types of sperm defects described for beef bulls 30 .The most serious morphological abnormalities were proximal droplets and loose heads.The high percentage of proximal droplets was 8.5% (Bull 446) but counts of 10-15% have been associated with decreased fertility 31 .Also, the highest percentage of detached heads was 18.9; although the bull can still be considered "fertile" with 30-40% of this defect, the ejaculate was considered infertile due to low viability (4%) (Bull 545).It has been described that the main causes of spermatozoa with detached heads are testicular hypoplasia, testicular degeneration and senescence of spermatozoa due to sexual inactivity.In this study the seminal samples of bulls were evaluated after the first ejaculation 30 .Our results were below threshold sperm values for compensable or uncompensable traits and provide evidence for the association between the sanitary status of seropositive-latently infected bull with BHV-1 and sperm morphology in Brahman bulls.
In general, it is accepted that bulls with values for progressive motility greater than 50% and normal morphology greater than 70% should be above the average fertility calculated for the sire group 29,30 .In this sense, the highest values for normal morphology, progressive motility and viability are positively associated with a higher non return rate than the average fertility in the group of bulls, but lower values of the seminal parameters expected as satisfactory for an ejaculate can be due to variability and do not indicate that the bull is subfertile.However, they do suggest that the bull with lower values should not be used for natural mating or as a semen donor in AI programs at the time 29,32 .Consequently, in this study high variability in sperm viability was observed between bulls with values that ranged from 4 to 93.8%.Particularly, in the group of bulls positive for BHV-1, it was observed that 50% (3/6) did not exceed 50% of viability, in accordance with what was reported by El-Mohamady et al. 33 although it indicated that BHV-5, a virus similar to BHV-1, does not induce functional and morphological damage of infected spermatozoa.In this study the type of morphological abnormalities was very diverse and the percentage of abnormalities decreased in the second sample (14.4%; range 3.8-23.4%)compared to the first evaluation (35%; range 22.6-51.6%),which suggests an effect of the sampling period rather than the effect of the virus on semen quality, contrary to what was reported by El-Mohamady et al. 7 because it was probably due to sperm turnover, in accordance to that indicated by the ABBA (2015).Recently, it has been reported in Brahman bulls that sexual rest decreases sperm quality 34 .More research is necessary on the issue of dead sperm in inactive/ rested bulls and exploring the possibility that sperm quality is better in the second sample.
In reference to the CASA variables, it seems that they do not contribute to differentiate seminal quality with a discriminant effect 35 .The differences observed in values of semen parameters for each bull, considering the sanitary status and the evaluations (0 and 30 days), show considerable variability and no clear associations were observed 2 .Furthermore, some authors suggested that BHV-1 does not affect semen quality 33,36 .Our findings indicate that although there were differences between sanitary status of Brahman bulls in kinematics variables, this did not appear to have a biological relationship with velocity, linearity, or undulation of sperm cells.
Based on CASA-Mot variables, our results identified four sperm subpopulations in Brahman bulls.The differences between sperm subpopulations (P < 0.05) suggest that sperm kinematics patterns with very active but non-progressive movement or rapid and progressive sperm movement may be indicative of early capacitation and hyperactivated-like motile spermatozoa.Other sperm subpopulations contained poorly motile sperm, and others presented non-progressive sperm and low velocity but greater progressive kinematics traits which highlights the heterogeneous nature of sperm ejaculates of fresh semen in Brahman bulls 18,37 .The analysis of sperm subpopulations by sanitary status provides evidence that the pattern of high velocity and progressive motility is associated with BLV − /BHV-1 + and BLV − /BHV-1 − , this suggestes that seropositive-infected bulls with BLV and BHV-1 should be discarded as sires.
The evaluation of seminal quality is very complex.Apparently, the variables that best discriminate between the quality of the ejaculates are normal morphology and sperm viability.The CASA-Mot variables do not reveal a clear trend with regards to the quality of semen samples.Taken together, the sanitary status of bulls did not show a clear change in the trend of the variables that define semen quality.However, the results indicated that the Table 3. Sanitary status and percentages of sperm anomalies in Brahman bulls (n = 10) at the end of the study (30 days).SS sanitary status: 1 = BLV + /BHV-1 + , 2 = BLV − /BHV-1 + , 3 = BLV − /BHV-1 − .www.nature.com/scientificreports/bulls infected with BLV or BHV-1 produce semen with variable quality, and that in most of the bulls the semen was free of the presence of the virus, but the use of bulls infected with BLV can generate semen contaminated in an intermittent way.The sperm subpopulation analysis indicated that seropositive-infected bulls with BLV and BHV-1 should be excluded from beef cattle farms.

Materials and methods
This study was performed following ethical principles and with the approval of the

Animal handling
During the experimental period, the animals were handled in an area of 18 hectares, a rotational system of seven apartments was established, with six days of occupation and 42 days of rest, with Ratana grass (Ischaemum indicum) and Mombaza (Panicum maximum).Animals were fed with grass without restrictions of access to forage consumption and they were supplied with salt, minerals and water ad libitum.The sanitary management of the bulls were an external dewormer (e.g., Doramectin), water-soluble vitamins (Catosal®) and minerals (Matsuda breeding Top Line).General inspection of the bulls was made, focusing on anatomical aspects such as conformation of legs and feet, body condition, coloration of mucous membranes and particularities of each animal.The condition of the external reproductive organs, such as the testicles, scrotum, foreskin, and penis were also examined.

Semen collection and processing
Semen was collected by electroejaculation using the Pulsator V® equipment (Lane Manufacturing Inc, Denver, CO, USA) as described by Montoya-Monsalve et al., and Víquez et al. 2,37 , under a collection program of paired semen samples.Semen was obtained twice, 30 days apart, early in the morning from ten bulls in sexual rest.
Semen was collected directly from the penis using a sterile graduate collection tube to assess the volume, each ejaculate was diluted at a ratio of 1:1 (vol/vol) using a commercial egg yolk-free diluent (OptiXcell 2-IMV Technologies, L' Aigle, France).Prior to semen dilution, a 1 mL aliquot of raw semen was placed in 1.5 mL sterile microtubes.Semen samples were assessed for gross motility, by placing 20 µL of fresh raw semen on a prewarmed slide at 37 °C, and for concentration, employing a bovine photometer AccuRead® (IMV Technologies, L' Aigle, France) at 530 nm wavelength.Diluted semen samples were placed on a heating plate at 37 °C and then transported in a polyurethane cooler at 25 °C.Subsequently, the raw semen samples were frozen at − 20 °C until nucleic acid extraction.Viability, percentage of sperm abnormalities and sperm morphology were assessed by preparing a smear of the semen sample using 10 µL of diluted semen stained with eosin-nigrosin as described by Kumar et al. 38 and at least 200 spermatozoa were counted to determine the percentage of abnormalities 22 .All samples were coded in such a way that the technician who performed the analysis did not know the number of the bull, the number of the ejaculate, or which ejaculate belonged to a particular bull.
After semen extraction, blood samples with and without anticoagulant were taken by coccygeal vein puncture using the Vacutainer system, and placed in a cooler at 4 °C.An aliquot of 1.5 mL of blood with anticoagulant was deposited in 1.5 mL microtubes and frozen at − 20 °C until the extraction of nucleic acids.Serum was separated by centrifugation at 3500×g 5 min, an aliquot of 1.5 mL was taken and stored in 1.5 mL microtubes at − 20 °C until assessment.

Assessment of sperm variables
Reusable counting chambers with cover slide for the analysis of motility and kinematics variables (Spermtrack® 20-micron depth, Proiser R + D, S.L., Paterna, Spain) were used after pre-warming them to 37 °C39 .After thorough mixing of the diluted semen samples, 2.7 µL of diluted semen were placed in the counting chamber tracks by drop displacement according to Valverde et al. 40 .Analyses were conducted with the CASA-Mot system ISAS®v1 (Integrated Semen Analysis System, Proiser R + D, Paterna, Spain) fitted with a video-camera (Proiser 782 M, Proiser R + D), a frame rate of 50 frames per second (fps) and a final resolution of 768 × 576 pixels as described by Valverde et al 41 .The camera was attached to a microscope UB203 (UOP/Proiser R + D) with a 1X eyepiece and a 10X negative-phase contrast objective (AN 0.25), and an integrated heated stage maintained at 37 ± 0.5 °C.

Computerized kinematic analysis
CASA analyses were performed recording seven microscope fields with a total of at least 600 cells per sample.The CASA-Mot variables assessed in this study included: straight-line velocity (VSL, µm•s −1 ), corresponding to the straight line from the beginning to the end of the track; curvilinear velocity (VCL, µm/s), measured over the actual point-to-point track followed by the cell; average path velocity (VAP, µm/s), the average velocity over the smoothed cell path; amplitude of lateral head displacement (ALH, µm), defined as the maximum of the measured width of the head oscillation as the sperm swims; beat-cross frequency (BCF, Hz), defined as the frequency with which the actual track crosses the smoothed track in either direction; motility (%), the percentage of total motile cells and progressive motility (%), corresponding to spermatozoa swimming rapidly forward in a straight line (assessed as straightness index ≥ 45%; VAP ≥ 25 µm/s).Three progression ratios, expressed as percentages, were calculated from the velocity measurements described above: linearity of forward progression (LIN = VSL/ VCL•100), straightness (STR = VSL/VAP•100), and wobble (WOB = VAP/VCL•100) 17 .

Serologic assessment using enzyme-linked immunosorbent assays
The sanitary status of the bulls for BLV and BHV-1 was determined in the Veterinary Virology Diagnostic and Research Unit, School of Veterinary Medicine, National University, Costa Rica using ID Screen® BLV gp51 Competition and ID Screen® IBR gB Competition (ID.vet,Grabels, France) according to the manufacturer's instructions, respectively.www.nature.com/scientificreports/Antigenic titration by polymerase chain reaction Genomic DNA extraction was performed from 300 µL of blood with EDTA and from semen diluted 1:5 (vol:vol) using Wizard® Genomic DNA Purification Kit (Promega Corporation, Madison, USA) according to the manufacturer's instructions for DNA extraction from samples of frozen blood.The DNA obtained was stored at -75 °C until use.BLV proviral DNA was identified by semi-nested PCR.In the first round, a set of primer targeting on gene encoding gp51 (envelope protein) were used, the forward primer was Env1: 5′-ATG CCY AAR RAA CGA CGR YCCCG-3′ (4826-4848), and the reverse primer was Env3: 5′-GAG GGC GAG RCC GGG TCC AGAGC-3′ (5550-5528), while the second round was performed with the forward primer Env2: 5′-CAA GGG CGR CGC CRR TTY GGAGC-3′ (5108-5124) and the reverse primer Env3 described above.The final volume of each PCR reaction was 20 µL and the mix consisted of 10 µL of Phusion U Multiplex PCR Master Mix #F-562S (Invitrogen, USA); 1.2 μL (0.3 μM/μL) of each primer pair Env1(+)/Env3(−) in the first round and Env2(+)/Env3(−) in the second round, 6.6 μL of RNase-free water (Invitrogen, USA) and 1 μL of proviral DNA in the first reaction or 1 μL of the initial PCR mixture of the first reaction after amplification in the second round to complete the volume for nested PCR.The amplification protocol was the same in both PCRs, as follows: initial denaturation at 98 °C for 3 min, 35 cycles of denaturation at 98 °C for 20 s, annealing at 70 °C for 30 s, and extension at 72 °C for 15 s, followed by a final extension at 72 °C during 5 min and maintenance at 4 °C thereafter.DNA from strain 280918-4418-CR isolated in Costa Rica was used as positive control (GenBank accession number: MN830810), while Rnase free water (Invitrogen, USA) was used as negative control.The expected products for BLV were 725 and 443 bp when the Env1(+)/Env3(−) primer pair was used in the first round or when the Env2(+)/Env3(−) primer pair was used in the second round, respectively, according to Khudhair et al 42 .
BHV-1 was identified using a set of primers that were based on region of the gE gen of BHV-1.The forward primer was BHV1-gE3-883F: 5′-CAA TAA CAG CGT AGA CCT GGTC-3′ (122,739-122,760), and the reverse primer was BHV1-gE3-989R: 5′-GCT GTA GTC CCA AGC TTC CAC-3′ (122,844-122,824) (Wernike et al., 2011).The final volume of each PCR reaction was 20 µL and consisted of 10 µL of Phusion U Multiplex PCR Master Mix #F562S (Invitrogen, USA), 1 µL of each primer IBR-IgE-F(+)/IBR-IgE(−) (0.25 µM/µL), 7 µL of Rnase-free H2O (Invitrogen, USA) and 1 µL of DNA.The amplification protocol consisted of an initial denaturation cycle at 98 °C for 2 min; 35 cycles of denaturation at 98 °C for 20 s, annealing at 63 °C for 30 s and amplification at 72 °C for 15 s, followed by a final amplification at 72 °C during 5 min and holding at 4 °C thereafter.DNA from the Colorado strain contained in the commercial vaccine EXPRESS® FP 10-HS (Boehringer Ingelheim Vetmedica, Germany) was used as a positive control.The expected product was 106 bp.All reactions were carried out in the 2720 thermocycler (Applied Biosystems, USA).The products amplified by PCR were stained with GelRed Nucleic Acid Stain 10000X in Water (Biotium, USA) and separated on a 2% agarose gel using the Fotodyne Model N°1-1430 electrophoresis equipment at 80 V during 40 min.Once separated they were visualized in the imager and ultraviolet illuminator (UVP BioDoc-1tTM Imaging System/FirstLight UV illuminator, USA).

Statistical analysis
A normal probability plot was used to assess normal distribution.The data obtained for the analysis of all sperm variables were assessed for homoscedasticity by using Levene test.Further, sperm variables were analyzed using the Generalized Linear Mixed Models (GLMM).The response variables were semen volume, total and progressive motility, swimming patterns (kinematics of spermatozoa), sperm concentration, normal and abnormal sperm (%), and viability (%).A normal distribution with an identity link function was assumed for all response variables.The quality of the seminal samples was categorized by analyzing the percentage of three variables (progressive motility, normal morphology, and viability when information was available), as follows: satisfactory seminal quality (S = total progressive motility > 50%, normal morphology > 70% and viability > 50%), doubtful (D = if at least one of the three variables did not reach satisfactory values) and unsatisfactory (U = if three or at least two of the variables did not reach satisfactory levels 19,20,29 .The association between semen quality with the sanitary status was evaluated descriptively.

Multivariate procedures
A subset of data was created with the means per ejaculate of all eight morphometric variables.Multivariate procedures were performed to identify sperm clusters from this subset of sperm morphometric data.All the values for the morphometric variables were standardized to avoid any scale effect.A principal factor analysis (PFA) was performed on these data to derive a small number of linear combinations that still retained as much information as possible from the original variables.Prior communalities for this analysis were estimated from the maximum absolute correlation coefficient between each variable and any other.The number of principal factors (PF) to be extracted was determined from the Kaiser criterion, namely by selecting only those with an eigenvalue > 1.The KMO (Kaiser-Meyer-Olkin) statistic was also obtained 43 as a measure of dataset adequacy for factor extraction.As a rotation method, the varimax method with Kaiser normalization was used 44 .
Further, an analysis was conducted to classify the ejaculates into a reduced number of clusters, based on scores obtained from factor analysis.This was accomplished in two-step phases, combining hierarchical and non-hierarchical clustering procedures with the sperm-derived indices obtained after the PFA.First, factor scores for all ejaculates were clustered hierarchically using the Ward Minimum Variance method 45 .From this analysis, an optimal number of clusters was determined based on criteria such as the Cubic Clustering Criterion (CCC), Pseudo-T, Pseudo-F, and partial R2.Second, the optimal number of clusters obtained in the previous analysis was used as the target number of clusters in a non-hierarchical K-means model cluster analysis 46 .
This classifies the spermatozoa of the data set into a small number of cluster (subpopulations) according to their head dimensions or kinematics variables, as has been described previously by Barquero et al. 47 .This analysis allowed the identification of sperm subpopulation and the detection of outliers.Analysis of variance (ANOVA) was further applied to evaluate statistical differences between clusters for all morphometric variables.Furthermore, an interaction effect between sanitary status and sperm subpopulation was considered.The threshold for significance was defined as P < 0.05.Multiple comparisons between cluster means were performed by the Tukey adjustment.Results were presented as mean ± standard error of the mean.All data were analyzed using the IBM SPSS statistical program, version 23.0 for Windows (SPSS Inc., Chicago, IL, USA).

Table 2 .
3anitary status, means of sperm parameters and semen quality of Brahman bulls (n = 10) at the end of the study(30 days).PM progressive motility, TM total motility, VCL curvilinear velocity, VSL straightline velocity, VAP average path velocity, ALH amplitude of lateral head displacement, LIN linearity of forward progression, STR straightness, WOB wobble, BCF beat-cross frequency, S satisfactory, D doubtful; U unsatisfactory. 1 Seropositive-infected bull with BLV and BHV-1.2Seropositive-latentlyinfectedbull with BHV-1.3Negativebull to BLV and BHV-1.

Table 4 .
Eigenvectors of principal components (PCs)* for the sperm kinematic for bull.VCL curvilinear velocity, VSL straight-line velocity.VAP average path velocity, LIN linearity of forward progression, STR straightness, WOB wobble, ALH amplitude of lateral head displacement.PC1 principal component designated "linearity", PC2 principal component designated "velocity".Var Exp variance explained in each PC (%).Total variance explained = 84.51%.*Expresses the more important variables in each PC.Only eigenvectors > 0.6 are presented.

Table 7 .
Percentage of each subpopulation (SP) in reference to the total number of spermatozoa the different sanitary status (SS).n number of spermatozoa.SS 1: BLV + /BHV-1 + ; SS 2: BLV − /BHV-1 + ; SS 3: BLV − /BHV-1 − a-d Different superscripts indicate significant differences; P < 0.05.Vol:.(1234567890)Scientific Reports | (2023) 13:18659 | https://doi.org/10.1038/s41598-023-45981-9 Committee of Centro de Investigación y Desarrollo de la Agricultura Sostenible para el Trópico Húmedo at the Costa Rica Institute of Technology (CIDASTH-ITCR) according to section 01/2019, article 1.0, DAGSC-074-2019.The care of the animals during the research complied with the animal welfare guidelines of Costa Rica Institute of Technology.Bulls were grazed and sampling procedures were reviewed and approved by the Production Committee of the Agricultural Production Program of the School of Agronomy at the Costa Rica Institute of Technology (PPASC-106-2021).A statement about plant ethics is not relevant because no experimentation is done on plants.The study was carried out in compliance with ARRIVE guidelines (https:// arriv eguid elines.org/; accessed on 8 November 2021).